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The data available for short-lived p-nuclides are used in an open nonlinear model of the 
chemical evolution of the Galaxy in order to discuss the origin of extinct radionuclides, the 
stellar sources of p-nuclides, and the chronology of solar system formation. It is concluded 
that the observed abundances of ®^Tc, ®®Tc, ®^Nb, and ^^®Sm in the early solar system are 
consistent with nucleosynthesis in type II supernovae during continuous chemical evolution 
of the Galaxy and a subsequent short isolation of the presolar molecular cloud from fresh 
nucleosynthetic inputs. However, further work on supernova models is needed before 
p-radionuclides will comprise reliable cosmochronometers. Despite these limitations, we 
argue that niobium-92 can be used to test whether the rp-process contributed to the 
synthesis of light p-nuclides in the Mo-Ru region. 

1. Introduction 

Extinct nuclides are radioactive nuclides that were present in the early solar system 
but have now decayed below detection levels. Their prior presence can be inferred from 
abundance variations of their daughter nuclides in solar system formation remnants. The 
origin of short-lived nuclides in the early solar system is the subject of much debate. 
They may have resulted from (z) inheritence from the chemical evolution of the Galaxy 
[ ^ H, {ii) injection from a nearby giant star that might have triggered the protosolar 
nebula into collapse [ , and {Hi) irradiation within the solar system by intense stellar 

flares [ I i i 0 i- In order to investigate the formation of the solar system, one can 
proceed either forwards or backwards in time. The abundances of short-lived p-nuclides 
in the interstellar medium at solar system birth can be derived from modeling of galactic 
chemical evolution [ 2, g] and stellar nucleosynthesis mm- The abundances of short-lived 
p-nuclides in the early solar system can also be derived from laboratory measurements of 
solar system formation remnants. Gomparison between the predicted and the observed 
abundances provides unequalled information on the origin of extinct radionuclides and 
solar system birth. 
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Table 1 

Abundances of short-lived p-nuclides. R is the initial ratio in the solar system derived from meteorite 
measurements, P is the production ratio in supernovae, 5 ?ess is the remainder ratio in the early solar 
system calculated as R/P, and SRism is the remainder ratio in the interstellar medium derived from GCE 
modeling (see text for details and references). 


Nuclide 


R 

P 

SRbss 

5Rism 

Technetium- 9 7 

y'Tc7«Ru 

< 4 X 10““ 

4.4 ± 1.3 X 10“^ 

< 1.3 X 10“^ 

1.2 ±0.3 X 10-® 

Technetium-98 

9®Tc/®®Ru 

< 8 X10“® 

7.3 ± 2.5 X 10“® 

< 1.7 X 10-2 

1.9 ±0.5 X 10-® 

Niobium-92 

92Nb/92Mo 

2.8 ±0.5 X 10“® 

1.5 ± 0.6 X 10-® 

1.9 ±0.8 X 10-2 

1.6 ±0.4 X 10-2 

Samarium-146 

i46Sm/i44sm 

7.6 ± 1.3 X 10-® 

1.8 ±0.6 X 10-1 

4.2 ± 1.6 X 10-2 

4.6 ±1.1 X 10-2 


2. Abundances in the early solar system 


Four short-lived p-nuclides might have been alive in the early solar system; which 
EC-decays to ®’^Mo (mean life of 3.8 Ma), ®®Tc which /^-decays to (6.1 Ma), ®^Nb 
which EC-decays to ®^Zr (50.1 Ma), and ^^®Sm which a-decays to ^'‘^Nd (149 Ma). Among 
these, only upper limits could be derived for ®^Tc [ 00 and ®®Tc [ 0 , while the 
abundances of ®^Nb [ [0 and [0 are hrmly established. The initial 


abundances normalized to a stable neighbor nuclide synthesized by the same process are 
compiled in Table 1. 


3. Stellar production ratios 

The origin of p-nuclides still poses some puzzles. The currently most favored production 
process is by photodisintegration of seed nuclides, this is called the 7 -process. However, 
as discussed below, there still remains the problem of Ending the proper site for the 7 - 
process. Furthermore, it is unknown whether the 7 -process is the only type of process 
that contributed to p-nuclide production. In the following we discuss two production sites, 
type II and type la supernovae. The feasibility of an additional scenario, the rp-process, 
will be discussed in Sec. |^. 


3.1. Type II supernovae 

Currently, core-collapse supernovae of type II (SNII) are strongly favored to be the site 
for the so-called p- or 7 -process outlined above. The necessary temperatures of about 
2 < Tg < 3 can be reached during hydrostatic and explosive neon/oxygen-burning [ 

|2D| , ^ pm , ^ |2^ , 0 . However, modern models still fail to reproduce the full range of 
p-nuclides although they succeed in producing substantial amounts in several mass ranges. 
More pronounced is the underproduction in the Mo-Ru region. For the initial production 
ratios ejected from the site of the p-process and used in our chemical evolution model 
described in Sec. 4, we have taken ratios from a recent study of nucleosynthesis in SNII 
[ [26[]. While p-nuclides with A > 100 — 110 are produced in amounts in rough agreement 
with solar abundances, for the Mo-Ru region the new model also exhibits underproduction 
features similar to those found in previous studies 100 . The reasons for this remain 
unclear. Details of the production are strongly dependent on stellar physics, such as 
the convection treatment used. Furthermore, the majority of the photodisintegration 
rates acting in this region are derived from theory. However, the experimental neutron 
capture data [ ^ as well as the rare data on proton capture in this region [ 0 agree well 
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with the predictions [ |2^, 0 . Additionally, an increased amonnt of seed nnclei wonld 
be needed in order to acconnt for the missing abnndances [ |^. For onr pnrposes, the 


absolnte prodnction factors are not important as we nse abnndance ratios. Obvionsly, in 
the following we have to assnme that those ratios will not change when the reason for the 
nnderprodnction is identihed. We nsed the SNII yields for three progenitor masses, 15, 21, 
and 25 Mq. These yields were integrated over the initial mass fnnction IMF, describing 
the birth rate 0 of stars with different initial mass m (</> oc dN/dm). For stars with initial 
mass higher than Mq, the initial mass fnnction follows a power law of the mass (j) oc m“", 
where a = 2.3 ±0.7 [1^. In order to calcnlate the IMF-averaged yields, the mass interval 


12.5-27.5 was further divided into three intervals 12.5-17.5, 17.5-22.5, and 22.5-27.5 Mq. 
The yields of the 15, 21, and 25 Mq mass stars were then weighted by the IMF integrated 
over the three mass intervals. The IMF-integrated SNII production ratios are compiled 
in Table 1. 

3.2. Type la supernovae 

Additional sites for the p-process have mainly drawn interest because of the problem 
with producing the Mo-Ru region. It has been proposed that the temperatures required 
for the photodisintegrations of the p-process could be reached in the C-rich zone of type 
la supernovae (SNIa) [ Early claims that the problem in the Mo-Ru region could also 
be alleviated have not been substantiated by making use of more realistic astrophysical 
models [ ^ |^. Howard and Meyer [ ^ End that they cannot produce p-nuclei in suffi¬ 
cient quantities to make a significant galactic contribution unless the s-process seeds are 
enhanced by several orders of magnitude. Furthermore, they observe that the produc¬ 
tion trends with respect to mass number are similar to those observed in SNII models. 
This can be understood by the fact that the burning during accretion on the C/0 core 
proceeds in a similar manner as the late burning stages of SNII progenitors. The actual 
photodisintegration process is the same. The above results seem to support the view that 
either SNIa do not contribute to the p-ratios or, if they do, the ratios are not signihcantly 
changed from those obtained with SNII models. 


4. Open nonlinear galactic chemical evolution 


The chemical evolution of the Galaxy in the solar neighborhood is governed by several 
mechanisms [ S H, 0. We assume that the surface density of the galactic disk in the 
solar neighborhood E grew with time from zero to its present-day value as a result of the 
accretion of low-metallicity gas. The infall rate is modelled as a normal density law Q 
with standard deviation equal to the mean r, where r is the infall time-scale [ ^ , 


dE 

dt 




( 1 ) 


The net evolution of the gas surface density Eg is a balance between star formation and 
disruption on one hand, and infall of low metallicity gas on the other hand. The rate of 
star formation depends on the gas density (Eg) according to Schmidt’s law. 


dt 


-uEl + ^gir). 


( 2 ) 
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For a nuclide for which the instantaneous recycling approximation is valid, its concentra¬ 
tion in the interstellar medium Z evolves through time as a result of stellar nucleosynthesis 
(y is the yield dehned as the ratio of newly created matter of that nucleus to newly cre¬ 
ated mass of permanent stellar remnants), dilution by infall of low-metallicity gas, and 
radioactive decay. 


dt 


yu.Ey‘-^|hl(Z-Z5)-AZ 

Zjg 


(3) 


Three parameters (.^, r, and a;) are enough to describe the chemical evolution of the 
Galaxy in the solar neighborhood and to characterize its most important features. In 
order to constrain these parameters, we proceeded like Sommer-Larsen [^. The present 


disk surface density, gas mass fraction, and metallicity in the solar neighborhood are 
known. Knowledge of these boundary conditions allows a straightforward derivation of 
the dependence of a;, and yo upon r. The G-dwarf metallicity distribution is used in 
order to estimate r and all other parameters which depend on it using a minimization 
algorithm. 

The exponent of the Schmidt power law is estimated to be 1.4 ± 0.15 


initial disk density is taken to be zero and its present-day value 45 ± 5 Mq pc 
The gas density is estimated to be 13 ± 4 Mq pc“^ 


g. The 


-2 


39 . The metallicity at solar system 


birth is [0/H] = 0. The metallicity of the infalling gas is assumed to have been constant 
through time, [0/H] = —1 [ ^ ^]. The age of the Galaxy is estimated to be 13.2 ± 1.5 
Ga [ The G-dwarf metallicity distribution in the solar neighborhood is taken from 
Haywood [ The model that best hts the G-dwarf metallicity distribution (Fig. 1) 
corresponds to = (56,4.8,0.065). 

Simulations provide a reliable means of estimating uncertainties on model parameters. 
Several sets of possible observations are hrst simulated. Then, each simulated set is 
used to compute a set of possible model parameters. Finally, this set of possible model 
parameters is examined for its distribution. The exponent of the Schmidt power law, the 
present surface density of the Galaxy in the solar neighborhood, and the gas density were 
simulated as normal density laws. Residuals between the best £t G-dwarf metallicity 
distribution and the observed distribution were calculated. These residuals were then 
resampled with repetition and added to the best £t distribution. This bootstrap procedure 
permitted the simulation of a set of possible G-dwarf metallicity distributions. Using the 
htted parameters, it is then possible to calculate the theoretical abundance of a p-nuclide 
in the interstellar medium at the time of solar system birth and its associated uncertainty. 


5. The remainder ratio 

The remainder ratio (3f?) is dehned as the ratio of the abundance of a radioactive 
nuclide to the abundance it would have if it were stable [ The remainder ratio in 
the interstellar medium at the time of isolation of the presolar molecular cloud from 
fresh nucleosynthetic inputs 3 ?ism can be evaluated within the framework of the galactic 
chemical evolution model that we have presented, 

3f?iSM= [ 

Jtq ^ Jtq ^ 


( 4 ) 
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Figure 1. G-dwarf metallicity distribution (density) in the solar neighborhood [Q. The [Fe/H] was 
converted to [0/H] using [0/H] = 0.70 ± 0.04[Fe/H] + 0.07 ± 0.02 [ Histogram: observations [ 
dotted line: closed-box model, solid line: best infall model. The parameters describing the best infall 
model are indicated at the top left corner. 


where 9 is the cycle number dehned as d9/dt = ^^(r)/Sg 
calculated in the previous section, it can be shown that 
remainder ratio takes the form, 


5B[]. Using the parameters 


for short-lived nuclides, the 


— «/(ATq), 


(5) 


with K = 2.7 ± 0.4, where the conhdence interval accounts for model uncertainties and 
Tq is the presolar age of the Galaxy 8.7 ± 1.5 Ga, /t=l corresponds to the closed-box 
model. Glayton [ estimated within a linear model that k must be within 2 and 4. The 
remainder ratio within the solar system at solar system birth (Ress can also be estimated 
using stellar production ratios and observed abundances in meteorites. 


»ESS = R/P, 


( 6 ) 


where R is the early solar system ratio of a radioactive nuclide to a stable nuclide produced 
by the same process and P is the production ratio. There might have been a delay between 
isolation of the presolar molecular cloud from fresh nucleosynthetic inputs and solar system 
formation, 

3f^ESS = ^^ISM C (7) 


where A is the decay constant and A is the free decay interval. If all short-lived p-nuclides 
were derived from the chemical evolution of the Galaxy, they should dehne a unique free 
decay interval. It is worthwhile to note that the deterministic chemical evolution model of 
the Galaxy probably breaks down for very short-lived nuclides because granularity might 
















6 


®'’Tc®®Tc “^Nb '“Sm 



Figure 2. The remainder ratio in the early solar system 3?ess is shown as a function of the remainder 
ratio in the interstellar medium at solar system birth 5 Rism for a variety of free decay intervals A. The 
points represent the positions of short-lived p-nuclides in this diagram (Table 1, see text for details). 


play a role, which would require a stochastic treatment [ The remainder ratios of 
®^Tc, ®^Tc, ®^Nb, and ^^®Sm are compiled in Table 1. As can be seen in Fig. 2, they are 
consistent with synthesis of short-lived p-nuclides in supernovae and derivation from the 
interstellar medium with a comparatively short free decay interval (A < 100 Ma). 


6. Origin of extinct radioactivities 

Assuming that there was no free decay interval between isolation of the presolar molec¬ 
ular cloud and condensation of the hrst solids in the solar system, our galactic chemical 
evolution model predicts that the initial ®^Tc/®®Ru and ®®Tc/®^Ru ratios at solar system 
birth were 5.3 ±2.0 x 10“® and 1.4 ±0.6 x 10“®, respectively. These values are consistent 
with the upper-limits derived from measurements (Sec. 1). 

There are a few very short-lived nuclides that cannot easily be explained by long-term 
chemical evolution of the Galaxy. These short-lived nuclides are ^Be, ^“^Be, ^®A1, ®^Mn, 
and ^^Ca. They may have been injected in the nascent solar system by a nearby giant 
star that might have triggered the protosolar cloud into collapse [ ^ . They may also have 
been synthesized within the solar system by local irradiation Both models 

face difficulties but the recent discovery of ^°Be [ in and possibly ^Be [ 0 in refractory 
inclusions favors the local irradiation model [ |^, H, 0, |^ . 

7. Niobium-92 and rp-nucleosynthesis in X-ray binaries 

As discussed previously, supernovae models fail to reproduce the absolute abundance of 
p-nuclides in the Mo-Ru region. This failure may be due to an inadequacy of the stellar 
models, involved nuclear physics, or might indicate that other stellar sources contributed 
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to the nucleosynthesis of p-nuclides. Interestingly, niobiuni-92 is a p-nuclide located in the 
mass range that shows a severe underproduction feature. Thus, this nuclide may provide 
useful constraints on p-nucleosynthesis. As indicated by other extinct r and p nuclides, 
it seems that the free decay interval was short, probably lower than a few hundreds of 
million years. Knowing the free decay interval, it is possible to do the reverse calculation 
and to compute the stellar ®^Nb/®^Mo production ratio using the observed abundance in 
the early solar system and the galactic chemical evolution model presented previously. We 
thus estimate that the stellar ®^Nb/®^Mo ratio was necessarily higher than 1.7±0.5 x 10“^ 
(A > 0) and was probably lower than 0.1 (A < 200 Ma). 

Another means of producing proton-rich nuclei is the rp-process 


46f| . In X-ray binaries. 


a thermonuclear runaway is ignited in the proton-rich layer accreted on the surface of a 
neutron star by mass flow from a companion star [ It has been shown that 

proton-rich nuclei can be synthesized in the subsequent explosive H and He burning. The 
rp-process path closely follows the proton drip-line. In the hnal freeze-out phase, the 
nuclides in the rp-path decay towards stability and form p-nuclides [ Recently, a 
dehnite endpoint of the rp-process was identihed which conhnes the possible production 
of p-nuclides to the mass region of about A < 110 [ Thus, this scenario would not 
interfere with the production in SNII but rather might account for the missing abundances 
in the Mo-Ru region. Among the isotopes considered here, the rp-process could contribute 
to the abundances of ®^Mo, ®^Tc, and ®^Ru. The nuclides ®^Nb and ®®Tc cannot be 
synthesized because they are shielded by stable nuclides. 

Niobium-92 is the only p-nuclide in the Mo-Ru region that cannot be synthetized in 
X-ray bursts but was alive in the early solar system. Thus, the abundance of this nuclide 
in the solar system may be used to test whether X-ray bursts signihcantly contributed to 
the synthesis of proton-rich nuclides. If X-ray bursts produced 90% of the p-nuclides in 
the Mo-Ru region and supernovae synthesized the rest, then the ®^Nb/®^Mo ratio expelled 
from stars to the interstellar medium would have been lowered by a factor of 10 relative 
to supernovae because X-ray bursts produce no ®^Nb. Thus, the apparent production 
ratio would have been 1.5 ± 1.0 x 10“^, a factor of 10 lower than the required minimum 
production ratio of 1.7 ± 0.5 x 10“^. This result speaks against a signihcant contribution 
of X-ray binaries to the synthesis of p-nuclides in the Mo-Ru region. Note that if X- 
ray binaries contributed to the synthesis of p-nuclides in the Mo-Ru region, they would 
overproduce nuclides in the A=104 region [^. The problem of the origin of the Mo and 
Ru p-nuclei therefore remains unsolved and clearly warrants further research. 
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